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Abstract. BV functions cannot be approximated well by piecewise constant functions, 
but we will show in this work that a good approximation is still possible with (countably) 
piecewise affinc functions. In particular, this approximation is area-strictly close to the 
original function and the -difference between the traces of the original and approxima- 
tion functions on the mesh can be made arbitrarily small. Necessarily, the mesh needs 
to be adapted to the singularities of the BV function to be approximated, and conse- 
quently, the proof is based on a blow-up argument together with explicit constructions of 
the mesh. We also establish a W^'^-error estimate for approximations by piecewise affine 
functions, which is used in the proof of BV-approximations and is also of independent 
interest. 
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1. Introduction 

Functions of bounded variation have important apphcations in many branches of math- 
ematical physics, among them optimization ^ free-discontinuity problems j3| (this and the 
previous reference also contains good introductions to the theory of BV functions) and hy- 



perbolic systems of conservation laws 10 . However, the theory of numerical approximations 
of such functions is not very well developed and indeed, because of the nature of singularities 
of BV functions, deeper analytical insight is requried. In this work we aim to contribute to 
such insight by considering the basic question whether BV function can be approximated 
well by (countably) piecewise affine functions. Before we start our investigation, though, we 
remark that any good approximation by piecewise-constant functions must fail, even though 
such ftmctions of course lie in BV, this will be shown in Proposition [3] below. 

For d > 1 an integer, let C M'* be a bounded Lipschitz domain. It is well-known 
that the space of piecewise afRne functions with countably many pieces is (norm-)dense in 
W^'P(ri; R™), where 1 < p < oo. Such a result is interesting for example in proving lower 
semicontinuity theorems, see for example [9 13 14 . Such an approximation result is also 



easily seen to be true in the space BV of functions of bounded variation, if we switch from 
norm convergence to the so-called strict convergence (see below). 

For some applications, however, one needs a better approximation result, in which also 
the trace differences over the boundaries of the mesh are controlled. This need becomes 
immediately obvious for example in the analysis of energy functionals in free-discontinuity 
problems, which typically involve a term measuring jumps over the discontinuity surface 
(see for example (Sl), and where, to produce good approximants, one needs to control the 
"jump error" . 
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Consequently, the purpose of this note is to prove the following approximation theorem, 



related to a weaker result in 17 and numerical and analytical studies |2l 5 6 21 , 25 



Theorem 1. Let ft C M.'^ be a bounded Lipschitz domain and u G BV(S1;IR'"). Then, for 
every £ > there exist a countable family TZ of (rotated) rectangles and simplices R d 
and V e W"'^'"'^(r2; M™) with the following properties: 

(i) The sets in TZ are non-overlapping and [C^ + \ IJT?.) = 0. 

(ii) For each R ^ TZ the restriction v\r is affine. 

(iii) \\u-v\\i^i^^,-i,^r^) + \{Du){n)~{Dv){n)\<e. 

(iv) I \u~v\dH''-^ <e. 

(v) w|ao = w|ao- 

Here, if = Cf" + is the Lebesgue-Radon-Nikodym decomposition of the measure 
/i, we define the measure (related to the area functional) 



(a.)(A):=|(m,£'^)|(A)^ y^^l+ ^"dx+|A*1(A) 

for every Borel set A C fi. 

A sequence (uj) C BV(f2;M") is said to converge strictly to u e BV(f2;]R'") if \uj - 
— ^ and \Duj\{Q) — )• |DM|(n) as j — )■ oo; it is said to converge (•)-strictly if even 
{Duj){fl) {Du){n). The latter convergence is the "natural" replacement for the strong 
convergence in BV(ri;M™) since smooth functions are (•)-strictly but not strongly dense in 
BV(ri;IR'"). In fact, we even have the following stronger result (see Lemma 1 of 15 or 
Lemma B.l in |7 for a proof under additional assumptions): 

Lemma 2. Let ft C M'' be a nonempty bounded and open set (no regularity assumption is 
imposed on dn) and let u G BV(f2;M™). There exists a sequence (vj) C (W^^^ nC°°)(17; E™) 
such that 

Vj — >■ u {•) -strictly and Vj ~ u e M™) for all j e N. 

If u £ W^'^{VL]M.™'), then we can in addition arrange that Vj — >■ u strongly in 

Many reasonable integral functionals with linear growth are continuous with respect to 
the (•)-strict convergence by Reshetnyak's continuity theorem [22! and extensions, see for ex- 
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ample Theorem 4 and the appendix of 16 ; lower semicontinuity in this situation is discussed 
in 11 . For some special examples of positively 1-homogeneous integrands see [20) . 

Thus, the main result of this paper establishes density of (countably) piecewise afSne 
functions in BV in the strict sense and such that additionally conditions (iv), (v) are satisfied. 
We remark that it is easy to satisfy conditions (i)-(iii) of Theorem [l] by simply mollifying u, 
which gives a strictly close smooth function (by Lemma [2| and then strongly approximating 
this smooth function by a piecewise affine function. Unfortunately, on this route one loses 
control of the sum of integrals in (iv), and it is precisely this condition, which is important 
for instance in applications to free-dicontinuity problems. 

In the course of the proof of the above theorem, we will rely on a result about approxi- 
mations in W^'^, which we factor out for its independent interest, see Theorem[7j It tells us 
that in W'^'"'^, we can prove actual approximation error bounds for any mesh, with the error 
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bound depending only the L^-modulus of continuity of the weak gradient of the approxi- 
mated function and regularity /uniformity constants of the mesh. Hence, in norm-compact 
subsets of W"'^'"'^ we get a uniform rate of convergence. 

Wc finally point out that the corresponding result for approximation by pure jump func- 
tions, that is, (countably) piecewise constant functions is not true — this also gives another 
reason why condition (iv) in Theorem [l] is of interest: 

Proposition 3. Let m e N 6e a natural number and assume that u € C^(B(0, 1), M™) 
is not identically zero. Then there can he no sequence of piecewise constant maps Uj € 
BV(B(0, such that Uj u {•)-strictly in BV. 

Proof. Let F: M™^'^ — > M be any continuous integrand satisfying F = on Vu(i?(0, 1)), 
which is a bounded set by assumption, and F{z) = \z\ for large values of \z\, hence F°° = 1. 
If Uj were piecewise constant BV functions such that Uj — >■ u (•)-strictly in BV, then by 



Reshetnyak's continuity theorem [22] (also see 24 and the appendix of [16]), we would have 

iDui I ^ \Du\ 



and, since Duj = D' 



where both convergences are weak* in CQ(i?(0, 1))*. But this is clearly impossible because 
Du^O. □ 

For the usual strict convergence, consider the continuous and positively 1-homogcncous 
integrand 

F:M2x2^]g^ |detA|i/2. 

Then, for any uj u strictly in BV(M^;M^), Reshetnyak's (original) continuity theorem 
implies 



dDui \ ,„ , * ^/ dDu \ ,„ , ^ 



However, if the Mj 's are pure jump functions, then rank | < 1 [-almost everywhere, 
but F is zero on the rank-one cone, and hence F{Du) = as a measure would follow, which 
clearly is not generally true. Thus it is not possible to approximate u with F(Du) 7^ by 
piecewise constant functions. 

This observation should be contrasted with the case of real-valued BV-functions (corre- 
sponding to m = 1) and the usual strict convergence, where a discretisation of the coarea 



formula allows approximation by piecewise constant functions, see 26 Theorem 3]. 

Finally, we note the following connection to quasiconvexity theory: Assume for a moment 
that we could find a piecewise constant sequence Uj u E BVo{il; M™) strictly (m > 2 and 
zero boundary values) and that there exists a positively 1-homogeneous rank-one convex, 
but not quasiconvex function F : M."^^'^ —>■ M. Then, by Reshetnyak's theorem and the fact 
that \D^Uj\-a..e. it holds that rank j'j^s",^. | — 1, 

\d\Du\J ' ' j^^J,, \d\D-Uj\J ' ^' 
> lim iD'ujUn) ■ F{0). 

But then F must be quasiconvex as well, a contradiction! 
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2. Approximation on a triangulation 

Let r be a rf-dimensional simplex in W^, that is, r is the convex hull of d + 1 affinely 
independent points vq, ■ ■ ■ , Vd of M.'^ (the vertices of r). Each point x G t admits a unique 
representation (the barycentric representation) as a convex combination of the vertices of r: 

d 

X = ^Xj{x)vj, 

3=0 

where Aj : t — t- [0, 1], j = 1, . . . , d, are affine functions with Yl'j=o = ^- If /• ^ 
is a continuous mapping, then the Lagrange interpolation of / on r is the (unique) afline 
mapping a: R'' — )• R™ that agrees with / at the vertices of r. Using the Aj's, it can be 
written as 

oo 
3=0 

Lemma 4. There exists a constant C = C{d,m), only depending on the dimensions, with 
the following property: Assume that f-.T^ is a & mapping with Lagrange interpolant 
a on T. Then, 

d 



IV/ - Va| dx < C(diamr)'^ ^ 



3=0' 

Proof. We start from the formula 

d 

^(a;-t>j)<8VAj(x) = -/ eM.'^'"^, 

3=0 



which is a consequence of VAj =0. Thus, 

/ d \ d 

Vf{x) - Va{x) = -V/(x) ( Y^ix - vj) ® VA,- J + ^(/(cc) - f{vj)) ® VA,- 

^ j=0 ^ 3=0 

j=0 ^-^0 ^ 



3=0 

Consequently, we find 



j \Vf{x) - Va{x)\ dx 

= Y.ffl + - • 1^ - ^^1' dH'^"'(^) ds dft, (2.1) 
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where we used (an elementary version of) the coarea formula for the last equality. Now, for 
fixed h G (0, 1) we change variables via y — $(a;, s) :— Vj + s{x — Vj), where (x, s) € {Xj = 
h} X (0, 1) and y £ t O {Xj > h}, for which we can estimate 

d'W'-Ux) ds = ^ jdi^mrr-^ 

^' |detD$| s'i-i/i(diamT) \y - Vj\d-^ ■ h 

since (diam r)s ^ |y — Furthermore, \x — Vj \ ^ (diamr)/! and hence 
IV-'fivj +s{x-Vj))\ ■ |x-ujf dn'^-^x) ds 

< C(diamT)2 f f 

Jo J{\j 

<C(diamr)^/ ^^^lif^ . h dy 
Plugging this in (2.1|, we conclude 

/|V/(a;)-Va(a;)|da;<C(diamr)^'V / / , ^^^'^^^A ■ h dy dh 



J{\j=h} 

1-1 r 

|VV(i;j + s{x - Vj))\ ■ dW^^^ix) ds 

,=h} 



This is the claim. □ 



Next we consider triangulations (we here always understand triangulations to be possibly 
infinite) . 

Definition 5. A triangulation o/il C M"^ is a family T = UjeN'^j 

(i) Tj — Gj+MjTQ for a fixed reference simplex tq C M'^ and some aj € M.'^ , Mj e GL{d). 

(ii) 17 = UjgN^j- 

(iii) Tj n Tfe = for all j ^ k. 

The triangulation T is called regular if there exists a constant a > with 
< I det Mj \<a for all j e N. 

We now have the following lemma: 

Lemma 6. Let u G BV(£';M™), where D C is a hounded polygonal domain. Further- 
more assume: 

(i) T — UjeN''j (possibly infinite) regular triangulation of D. 

(ii) T satisfies a hounded overlap property, that is, there exists a constant iiT G N such 
that every x € D lies in at most K different sets r + i?(0, diamr), where r G T. 

Then, there exists a piecewise affine function a: D ^ M™ with respect to the triangulation 
T and a constant C — C{d, m) such that 

[ \a\dx<C [ \u\dx and [ \\7a\ dx < C\Du\{D). (2.2) 
Jd Jd Jd 

Proof. We only show the second inequality, the first is similar. We also implicitly consider 
u to be extended by a constant to m G BY{D + B(0, i?); M") such that \Dii\{D + B{0, R)) < 
C\Du\{D), where R = max-,- diamrj-; this is possible for example by Lemma 9.5.5 in 18 
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Apply Lemma|4]to the regularized mapping / = ^Pg-ku, where v? S C^(M'^) with supp tp C 
i?(0, 1) and {'p)e{x) :~ e^^Lpixje) for e > to be chosen later. For integers 1 < fc, Z < d we 
have 

dki! = -{diip)E * {dku - z) for any z e M™. 

e 

Let r be a fixed simplex in the triangulation T. Consequently, we may estimate using 
Lemma [4j 



J |V/(a;) - \7a{x)\ dx 



S=0 



< ^ t(diam.)^s-||V^||,. / / \'^<x-y) z\ 

< ^l|V¥'llL-(diamr)^^ / \\/u{w) - z\ [ 



J(0,e) JBiO.e) \w + IJ - Vj 



aw. 



Id-l 



Here we remarkt that all constants can be chosen independently of t, because the triangu- 
lation is regular. Taking 

e = diamr, 

we note 

dy r dy f dx 



and so, 



/ |V/- Va| dx < Cinf /" |Vw-z|da;, (2.3) 

Jt ^ v/t+S(0, diamr) 

hence, setting z = and summing these integrals over all simplices t using the bounded 
overlap property of the triangulation, 



/ |Va| dx<C [ |Vu| + |V/| dx, 

Jd J D+B(0.dia.mT) 



ID J D+B{0,dia.mT) 

which immediately implies the claim. □ 

Finally, we remark that for W^'^-functions u, we also have convergence to the original 
function on uniform regular meshes and we can even quantify the speed of convergence in 
terms of the L^-modulus of continuity for Vu, defined as 

uj{h) :— sup / \u{x + w) — u{x)\ dx . 

\w\<hJD 

The result then is: 

Theorem 7. In the situation of Lemma^ assume additionally that u G W"'^''^(ri; M"*) and 
that the triangulation T is uniform, that is, there exists a constant 7 G (0, 1) such that 

"fk < diam t < k for all r G T, 

where 

k := sup diamr. 
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Then, there exists a constant C , only depending on the dimensions and the parameter of 
regularity of the mesh T , such that 



/C 
|Vw- Va| dx < ^uj{ik). 



It is well-known that ijj{h) — >■ as /i ^> (proof by approximation with smooth functions), 
but of course the speed of convergence is highly dependent on u, so in general there is no 
"rate" of convergence. However, by the Frechet-Kolmogorov characterization of compact 



sets in W ' , see for example 12 , there is a uniform rate of convergence for all functions u 



in a fixed compact subset of W ' (Jl; 



Proof. We abbreviate t+ r + i?(0, diamr) and plug z :— (Vm)t-+ :— into (2.3). Then 
we estimate 



j I V/ -Va\dx<C J I Vu - (Vm)^+ | dx 

<c/ / \\/u{x) ~\/u{y)\dy dx 



J T 



<Cf / \Vu{x)-Vu{x + w)\dwdx 

r+ "'B(0,diamT+) 



<Cf / \Vu{x)~Vu{x + w)\dxdw. 

< , , / / \V u{x) — V u{x + w)\ dx dw . 

7 JS(0,3fc) Jt+ 

Thus, summing over all r G T and using the bounded overlap property of the system 
{^*}reT. 



(0,3/c) JD 



|Vu(a;) — \7u{x + w)\ dx dw 



Using the fact that ||Vu — f\\'L^{D;W") can be made arbitrarily small, this finishes the 
proof. □ 

From the preceding theorem we immediately get convergence with a rate if Vu is more 
regular: For example, if w € W^^^(£';M™), then we have convergence with a linear rate (on 
regular uniform meshes). 

3. Proof of the main result 
We first establish the following local lemma: 

Lemma 8. Let Q he an open (possibly rotated) cube in , let Qq CC Q be a concentric 
similar open subcuhe, and u G BV((5;M'"). We further assume: 

(i) The inner cube Qq is the disjoint union of a finite number of similar (closed) rect- 
angles Si, ... , Sm that are translations of one another along a single axis parallel to 
one of the sides of Qq, and with the length rj of the short side of the Sj , it holds that 
dist(9(5, Qo) = rj (like in Figure^. 
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- (1/2 + 1/4 + 1/8)77 

- (1/2 + 1/4)77 




Figure 1. Triangulation used in the proof of Lemmajs] (77 = dist(Q'^, Qq)). 

(ii) We are given a function w G {Qo;M."^) that is affine when restricted to any 

Sj, and satisfies 

n „ 

J2 dH''-i < e\Du\{Q). 

Then, on A := Q \ Qq there exists a countably piecewise affine mapping v S W^'^(^;M™) 
satisfying v = u on dQ, v — w on OQq, and 

||t'||Li(A;R-) < C||ii||Li(A;R") and (Dv) {A) < C {C^ + {Du)){A) . (3.1) 
Moreover, if T is the triangulation of A corresponding to v, then 

[ \u-v\dn'^-^ <C\Du\iA). (3.2) 
Here, C = C{d,m) is a dimensional constant. 



Proof of Lemma^ We start by constructing a triangulation T = UjeN '^^ ^ Figurejlj 
This triangulation is of Whitney- type towards the outer boundary dQ, 

diamr ^ dist(T, dQ). 

Our T also has the property that the simplices match the elements Si, . . . , Sn at the inner 
boundary OQq, meaning that for each simplex t € T we have r H dQo C dSk for some k, 

At this stage we invoke Lemma[6j This allows us to find v G W^'^(^; M™) with 1; = u on 
dQ, V = w on dQo and such that p.l[) holds. For the latter note that |-\/l — \A\'^ — \A\\ < 1, 



hence the second statement in (3.11 follows from (2.2) 
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Now for r > and t & T we let t'^^^ denote the simplex obtained by dilating r by r about 
its centre of mass. In view of our construction of T, we can find an absolute r = r{c, d) > 1, 
such that for each r e T, a substantial fraction of r'^''^ touches the outside of Q; more 
precisely, we require that £''(t*^''' \ Q) > 26C'^{t) for all t e T and for some constant (5 > 0. 
We define for each t £T, t* ^ t^''^ \ Qq. Then, 

C'i{T* \ Q) = /:'^(t('-) \Q)> 25C'^{t) > SC^iT*) for ah r e T. 

For checking ( |3.2[ ) we note that we can extend u — v to M'' by zero, and that D{u — u) = 
Dv — Du for this extension since v ^ u on DQ. Since m — w = on r* \ Q we have by standard 
Sobolev and Poincare inequalities (see for instance Theorem 4.4.2 in [2^, the specific form 
of the constant can be derived by a simple scaling argument as in Remark 3.50 in [s]) 



\u - v\ AW^' < C (dianir)"' / \u~v\Ax + \D{u ~ v)\{t) 

< C ^(diamr)-i j \u-v\Ax + (\Du\ + \Dv\){t* n Q) 

<C{^^ + l){\Du\ + \Dv\){T*f^Q) 
<C{\Du\ + \Dv\){t* r\Q). 



Consequently, we get using ( |3.1| , the bounded overlap property of the system t* , t £ T, 
and the result of Lemma |6] (in particular the fact that its constant is universal) that 



f lu-vldH''^^ <cY,{\Du\ + \Dv\){t* DQ) 



<C{\Du\ + \Dv\)iA) 
< C\Du\{A). 

Of course, the constant C depends on the r we chose, but since the geometry of the trian- 
gulation is fixed, this can be chosen absolutely. □ 

We can now show our main result. 

Proof of Theorem^ We introduce the following notation: Let C = (—1, 1)"^ and for xq G M.'^, 
r > set Q{xo, r) :— xq + rC. For every unit vector n E S''^^ select a rotation P e SO{d) 
with Pei — n, where ei = (1,0,..., 0)^ e R''. Let Qn{xo,r) :— xq + rP(C), so that 
Qnixa,r) is an open cube with center xq, sidelength 2r, and two faces orthogonal to n. 
Step 1. For £'*-almost every xq E^ and < r < (i\st{xQ^d^) / \/d put 

u{xo + ry) - u{xo) 
r 

where u{x()) is the value of the precise representative of u at a:o (for this and other properties 
of BV- functions see [s]). We recall that for £''-almost all xq € Vl, u is approximately 
differentiable at xq and so, as r | 0, 

Ur->uo (•)-strictly in BV(C;M™) with uo{x) ■.^Vu{xf^)x, x E C, (3.3) 

with the approximate gradient Vm(xo) of u at Xq. The fact that one can indeed choose a 
sequence r„ J, such that — > uo (•)-strictly follows for example from Lemma 3.1 in [23| 
about the construction of strictly converging blow-ups, applied to the measure {Du^C^). 
Since the trace operator is strictly continuous, we also have that 

/ \ur - uo\ dW^-^ ^ Q as r i 0; (3.4) 
Jac 
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see for instance 16 , pp. 53-54, and the references given there. Let us denote the set of such 
points xo S by Gi. 
Now, for 

w{x) = u{xo) + '^u{xo){x — xo), xeQ{xQ,r), (3.5) 



from (3.3), ( |3.4[ ) we get by a change of variables that there exists < r(a;o) < 1 such that 
for r < r{xo) we have 

\u — w I dx < erC'^{Q{xo,r)) 



Q(xo,r) 

{Du)iQ{xo,r)) - {Dw)iQixo..r))\ < eC\Q{xo,r)), 

\u-w\A'H'^-^ <eC'^{Q[xo,r)). (3.6) 

dQ{xo,r) 

Step 2. Next, for |£)*M|-almost all xq g f], henceforth fixed, we have Ur := ^'""^^'^^d''^"'''"''' 
00 as r I 0, and, defining 

u{xQ + ry) - u^^^r ^ ^ 

rar 



where Uxo^r = j-Q (^^o r) ^ Lemma 3.1 of 23 we can find a sequence of r's going to 

(not specifically labelled) such that 

Ur uq strictly in BV(Q„(0, 2); M™) (3.7) 

and 

uo{y)^bij{yn), yeQ„(0,2), (3.8) 

where ^/^ : (— 1, 1) — )■ M is increasing and bounded. The fact that the blow-up limit uq can 



indeed be chosen one-directional can be proved using the rigidity result in Lemma 3.2 of 23 



(also see Remark 5.8 in loc. cit.), or using Alberti's rank one theorem |T], see Theorem 3.95 
in [3]. 

It is not hard to see that we additionally may assume that the sequence of r's is chosen 
such that \Du\{dQn{xQ, r)) — 0; we again refer to 16 , pp. 54-55, for these assertions. Denote 
the set of such points € by G2, and observe that {C''- + \Du\){il \ (Gi U G2)) — 0. 

Let N satisfy 

iV>"----^^^"^'(-^-^)^^>'^^(2)-^(-2)). 



For the function tp appearing in (3.8) we claim that we can require that the equidistant 
partition 

2 

'n' 



-\ = to < ti < ... < ti^ ~ 1 with tj+i — tj 



consists only of continuity points of i/i. This can be achieved as follows: Select < 6* < 1 
such that for the modified function i)g{t) := ^(t + 9) all the tj are continuity points. Since 
has only countably many discontinuity points, such 9 always exists. This corresponds to 
a blow-up sequence of the form 

u{xQ + r9n'Vry)-Uxo,T ^ n 1 ^ 

rar 

and in the following we need to replace Qn{.XQ,r) by Qn{xQ + r9n,r). We note that N still 
satisfies 
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In the following we will however suppress such a possible shift for ease of notation. 

Define </? as the piecewise affine function satisfying ip{tj) = 4^(tj) for each j = 0, 
and note that with (the rotation P as above), 

S, :=P[(i„i,+i) X (-1,1)^-'], i = 0,...,7V-l, 

we have 



\uo{y) ~ b^{y ■n)\dx^ 2''-'\b\ j ^\^-^\ds 



(3.9) 



3=0 



and 



J2 I \u^{y)-b^{yn)\dH''-\y) 

= l^lE / \Hyn)~v{yn)\m''-\y) 



< 



=0 -^dSj 



N 



J2 ^{tj+i) ' i^itj) 



< 



(3.10) 



Next, in view of (3.7 1 and our choice of partition points tj, we infer from the trace theorem 
that 

AT-l 



E 



\Ur — Uo\ dH 



d-1 



as r I 0. 



''9Sj 



(3.11) 



For a point xq G G2, the mapping x 1— > bip{x ■ n) defined above is piecewise affine. 
Split the the rotated cube Qn{xQ,r) into k rectangles Sj{xo,r) := xq + rSj and define the 
corresponding piecewise affine map 



'W{x) := Uxo,r + rurb ip 



X — Xq 



X e QnixQ,r). 



(3.12) 



Hence, changing variables in (3.7), (3.11) and using (3.9), (3.10) as well as the estimate 
I Y^l — |Ap — 1^11 < 1 it follows that there exists r{xo) > such that for r < r{xo), 

\u - It; I da; < er\Du\{QnixQ, r)). 



|(Du)(Q„(xo,r)) - (Dw)(Q„(xo,r))| < e{Du){Qn{xo,r)) + 2C{Qn{xo,r)), 

N-1 „ 

E/ \u-w\dn''-^ <e\Du\{Qn{xo,r)). (3.13) 

Step 3. For every G G := Gi U G2 we have so far constructed a (rotated) cube 
Qq — Qn{xo,r) with n and r depending on xq, and for all xq G G2 this cube is further 
subdivided into rectangles Sj{xo, r) [j = 0, . . . , N — 1). Now, for every such Qo we choose 
a slightly larger concentric similar cube Q = Q{xQ,r) Qq with the properties 

II"IIli(q\q^;R"^) < £||"IIli(Qo;R'")- + {Du)){Q \ Qo) < e{Du){Qo) (3.14) 

and in the case a;o G G2, we further require 

2r 



dist{dQ,Qo) = rj 



iV' 
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where rj is the length of the short side of the Sj {xq, r')'s. This can be achieved by subdividing 
the rectangles and correspondingly shrinking Q. 

Now we invoke Lemma [s] with Qo and Q and with w as in ( |3.5[ ) or (3.12) for xq G Gi or 
xo G G2, respectively. In particular, this yields a (countably) piecewise affine function vq 
with the properties stated in that lemma. 

Step 4- We next apply the Morse Covering Theorem [l9j (or see Theorem 5.51 in [s]) 
to cover {C^ + |I?u|)-almost all of il with (rotated) cubes Q from the above family. After 
subdividing the cubes Qo = Qo{xo,r) corresponding to points a;o e G2 into the rectangles 
Sj{xo:r) and the set Q{xQ,r) \ QoixQ,r) into the simplices constructed in Lemma [sj we 
thereby find a family TZ satisfying (i) in the statement of Theorem [T] 

For the remaining properties (ii)-(v), we write 

where an denotes the afhne map corresponding to the rectangle or simplex R d TZ (in par- 
ticular Vq = X]_RcQ ^-R ^'^^ froni the previous step). Because the rectangles/simplices 
are disjoint, the map v is well-defined and it is clear that (ii) is satisfied. For the remaining 
assertions we consider for j G N the mapping 

/ 



with 



u ^ 

\mR)>] ) 



Since the above sum is finite, we infer that Vj G BV(il;M''") and 



J2 {yaRC^^R + an^vRH'^-^^dR) 



c^(R)>^ 



Du^H.j+ {VaRC'^^R+{aR-u)®VRn'^-^^dR) 



C-^iR)>i 



Here, and i>r are the unit inner normals to Hj and R, respectively. Employing (3.6 1, (3.13 ), 
\Dv,\{n) ^ \Du\{Hj) + V (\^aR\C''{R)+ [ \aR ~ u\ dV."- 

1 V JdR 



£<*(i?)>i 



= \Du\{H,)+ J2 \yaR\C''iR) + 0{e){C''+\Du\)in). 

C''{R)>j 

Since \Du\{Hj) ^ as j ^ 00, we see that v G BV(f^;M™). 



Concerning (ui), we estimate using ( |3.6[ ), (3.13), (3.1), (3.14), 

\{Du)in) - {Dv)in)\ < Y,[\{Du){Qo) - {Dv){Qo)\ + {{Du) + {Dv)){Q\Qo) 

Q 

Q 

<eCY,{C'' + {Du)){Qa)+C\Z) 
Q 

< eCiC" + {Du)){ir) + C\Z), 
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where the summation is over ah cubes used to cover il; inside the sum Qq = Qn{xo, r) refers 
to the inner cube. The term C^iZ) originates from the additional Lebeseue measure for the 
cubes corresponding to singular points x E G2- Because they only occur on a £'*-negligible 
set, we can make this term disappear in the limit e — > 0. 

By a similar calculation, also (iv) holds (this time we use (3.2) instead of (3.1 1): 

\u-v\ dW^-^ 



E 

Ben 



E 



N-l 

/ \u-v\dn''-^ + Y, / \u-v\dn''-^ 

j=Q ''dSj{xo,r) ^^^Jot 



< 



+ \Du\)iQ) + C\Du\iQ^ixo, r) \ Qoixo, r)) 



<e{l + C){C'^ + \Du\){n). 

For (v) we only need to observe that Vj — u on dfl for all j E N and Vj —>■ u strictly, hence 
the trace is preserved. 

The only remaining part to check is whether the constructed mapping v is of class 
W^'^(r2; M™). The features we shall use here are that for every cube Q as before, v\q E 
W^'^((3; M™), t; = M on dQ, and that we may assume that \Du\{dQ) = 0. We have 

Dv ^ Lg + v\aQ vqW^-^^OQ). 

Q 

By assumption v\qq — u\qq and the latter coincides also with the outer trace of u on dQ 
since \Du\{dQ) = 0. Keeping in mind that X^qgq -"-Q ~ -"-f^ /I'^-almost everywhere, and 
hence in the sense of L^(M''), we find 

QeQ QeQ 



This concludes the proof. 



□ 
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